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Abstract 

This paper considers exponential utility indifference pricing for a mul- 
tidimensional non-traded assets model and provides two approximations 
for the utility indifference price: a linear approximation by Picard itera- 
tion and a semigroup approximation by splitting techniques. The key tool 
is the probabilistic representation for the utility indifference price by the 
solution of fully coupled linear forward-backward stochastic differential 
equations. We apply our methodology to study the counterparty risk of 
derivatives in incomplete markets. 
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1 Introduction 



The purpose of this paper is to consider exponential utility indifference pric- 
ing in a multidimensional non-traded assets setting, which is motivated by our 
study of counterparty risk of derivatives in incomplete markets. Our interest is 
in pricing and hedging derivatives written on assets which are not traded. The 
market is incomplete as the risks arising from having exposure to non-traded 
assets cannot be fully hedged. We take a utility indifference approach whereby 
the utility indifference price for derivative is the cash amount the investor is 
willing to pay such that she is no worse off in expected utility terms than she 
would have been without the derivative. 

There has been considerable research in the area of exponential utility indif- 
ference valuation, but despite the interest in this pricing and hedging approach, 
there have been relatively few explicit formulas derived. The well known one 
dimensional non-traded assets model is an exception and in a Markovian frame- 
work with a derivative written on a single non-traded asset, and partial hedging 
in a financial asset, Henderson and Hobson |21j, Henderson [T8j, and Musiela 
and Zariphopoulou [44] use the Cole-Hopf transformation (or distortion power) 
to linearize the non-linear PDE for the value function. This trick results in an 
explicit formula for the exponential utility indifference price. Subsequent gen- 
eralizations of the model from Tehranchi , Frei and Schweizer [TB] and [T7] 
have shown the exponential utility indifference value can still be written in a 
closed-form expression similar to that known for the Brownian setting, although 
the structure of the formula can be much less explicit. 

As soon as one of the assumptions made in the one dimensional non-traded 
asset breaks down, explicit formulas are no longer available. For example, if the 
option payoff depends also on the traded asset. Sircar and Zariphopoulou |47j 
develop bounds and asymptotic expansions for the exponential utility indiffer- 
ence price. In an energy context, we may be interested in partially observed 
models and need filtering techniques to numerically compute expectations (see 
Carmona and Ludkovski [II and Chapter 7 of [TIT). If the utility function is not 
exponential, Henderson 18^ and Kramkov and Sirbu |32j developed expansions 
in small quantity for utility indifference prices under power utility. 

In this paper, we study the exponential utility indifference price in a multidi- 
mensional setting with the aim of developing a pricing methodology. The main 
economic motivation for us to develop the multidimensional framework is to con- 
sider the counterparty default risk of options traded in over the counter (OTC) 
markets, often called vulnerable options. The recent credit crisis has brought to 
the forefront the importance of counterparty default risk as there were numer- 
ous high profile defaults leading to counterparty losses. In response, there have 
been many recent studies (for example, Bielecki, Crepey, Jeanblanc and Zargari 
[B] and Brigo and Chourdakis [5]) addressing in particular the counterparty risk 
of CDS. In contrast, there is relatively little recent work on counterparty risk 
for other derivatives, despite OTC options being a sizable fraction of the OTC 
derivatives market^ The option holder faces both price risk arising from the 
fluctuation of the assets underlying her option and counterparty default risk 
that the option writer does not honor her obligations. Default occurs when the 

^In fact, OTC options comprised about 10% of the $600 trillion (in terms of notional 
amounts) OTC derivatives market at the end of June 20101 whilst the CDS market was 
about half as large at around $30 trillion. 
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assets of the counterparty are below its liabilities at maturity (following the 
structural approach of Merton). In our setting, the assets of the counterparty 
and the asset underlying the option are non-traded and thus a multidimensional 
non-traded assets model naturally arises. A second potential area of application 
is to residual or basis risks arising when the asset (s) used for hedging differ from 
the assets underlying the contract in question (see Davis |13p. 

Our first contribution is to use the solution of fully coupled linear forward- 
backward stochastic differential equation to give a probabilistic representation 
for the exponential utility indifference price. Since the associated equations are 
linear (but with a loop), we call it a pseudo linear pricing rule. It is well known 
that the utility indifference price can be written as a nonlinear expectation of 
payoff under the original physical measure, and the nonlinear expectation is 
often specified by a backward stochastic differential equation with quadratic 
growth {quadratic BSDE for short). Several authors derive quadratic BSDE 
representations of exponential utility indifference values in models of varying 
generality - see Mania and Schweizer [S^, Ankirchner et al [1], Becherer [5], 
and Frei and Schweizer [TB] |T7] among others. In contrast to representing 
the utility indifference price by nonlinear expectation, we represent it as the 
linear expectation of payoff under some equivalent pricing measure. The pricing 
measure we choose may depend on the utility indifference price itself, so the 
pricing mechanism is a loop. What makes this transfer possible is the nonlinear 
Girsanov's transformation in the sense that the change of probability measure 
involves the solution of the equation, and the BMO martingale property of 
the solution, both of which essentially follow from the comparison principle of 
quadratic BSDEs (see for example Hu et al [53] )• This kind of pseudo linear 
pricing rule appeared in Proposition II of Mania and Schweizer [39], where they 
modeled the payoff as a general random variable. In contrast, as we specify the 
dynamics of the underlying assets and the payoff structure, a fully coupled 
linear FBSDE appears naturally. Since the associated equations are linear (but 
with a loop), Picard iteration can be employed to approximate the solution 
and the corresponding utility indifference price. This can be regarded as a first 
application of the pseudo linear pricing rule. 

Our second contribution is to develop a semigroup approximation for the 
pricing PDE by the splitting method, where we specialize to a Markovian set- 
ting and derive the semilinear pricing PDE with quadratic gradients satisfied 
by the utility indifference price. In our multidimensional setting, the Cole-Hopf 
transformation (as in the one dimensional model) can not be applied directly 
since the coefficients of the quadratic gradient terms do not match. Motivated 
by the idea of the splitting method ( or fractional step, prediction and correc- 
tion) in numerical analysis, we split the pricing equation into two semilinear 
PDEs with quadratic gradients such that the Cole-Hopf transformation can be 
applied to linearize both of them. Splitting methods have been used to con- 
struct numerical schemes for PDEs arising in mathematical finance - see the 
review of Barles [5], and Tourin ^50j with the references therein. Recently, 
Nadtochiy and Zariphopoulou [IS] applied splitting to the marginal Hamilton- 
Jacobi-Bellman (HJB) equation arising from optimal investment in a two-factor 
stochastic volatility model with general utility functions. They show their 
scheme converges to the unique viscosity solution of the limiting equation. 

The idea of splitting in our setting is as follows. The local time derivative of 
the pricing PDE depends on the sum of semigroup operators corresponding to 
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the different factors. These semigroups usually are of different nature. For each 
sub-problem corresponding to each semigroup there might be an effective way 
providing solutions. For the sum of these semigroups, however, we cannot find 
an accurate method. The application of splitting method means that instead of 
the sum, we treat the semigroup operators separately. We prove that when the 
the mesh of the time partition goes to zero, the approximated price will con- 
verge to the utility indifference price, relying on the monotone scheme argument 
by Barles and Souganidis any monotone, stable and consistent numerical 
scheme converges (to the correct solution) provided there exists a comparison 
principle for the limiting equations. The difficult part of applying the monotone 
scheme to our problem is the verification of consistency. This is overcome by 
the pseudo linear pricing rule for the utility indifference price. In contrast to 
the nonlinear expectation, where the Dominated Convergence Theorem may not 
hold, the representation is linear, so Dominated Convergence Theorem can be 
employed to verify the commute of limiting processes. This can be regarded as 
a second application of the pseudo linear pricing rule. 

Our third contribution is to apply the splitting method to compute prices of 
derivatives on a non-traded asset and where the derivative holder is subject to 
non-traded counterparty default risk. In contrast to the complete market Black 
Scholes style formulas obtained by Johnson and Stulz [5H] , Klein [53] and Klein 
and Inglis {30| , we show the significant impact that non-tradeable risks have on 
the valuation of vulnerable options and the role played by partial hedging. 

The paper is organized as follows: In Section [21 we propose our multidimen- 
sional non-traded assets model, and present our probabilistic representation for 
the utility indifference price. Based on such representation, we give the Picard 
approximation for the price. In Section [3l we specialize to a Markovian setting, 
and apply the splitting method to give a semigroup approximation for the util- 
ity indifference pricing PDF. We further apply the scheme to the counterparty 
risk of derivatives in incomplete markets in Section [Jj 

2 Pseudo Linear Pricing Rule for Utility Indif- 
ference Valuation 

Let yV = {W-^, • • • , W^) be a d-dimensional Brownian motion on a filtered 
probability space J^, {-T-j}, P) satisfying the usual conditions, where Tt is the 
augmented tr-algebra generated by {Wu : < u < t). The market consists of a 
traded financial index P, whose price process is given by 

—L=f,[dt+{a[,dWt}, (2.1) 

and a set of observable but non-traded assets S = {S^, ■ ■ ■ ,5'"), whose price 
processes are given by 

-± = ^^ldt+{al,dWt) (2.2) 

for i = 1, ■ ■ ■ ,n. (•, •) denotes the inner product in M.'^ with its Fuclidean norm 
II • ||. We have ^^f , ^^l e M, af = (af i, • • • , G and aj = (aji, • • • , al^) S 
K*^. There is also a risk- free bond or bank account with price i?t = 1 for < > 0. 
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This is equivalent to working in terms of discounted units and is without loss of 
generality. 

Our interest will be in pricing and hedging (path-dependent) contingent 
claims written on the non-traded assets S. Specifically, we are concerned with 
contracts with payoff at maturity T of g{S.) which may depend on the whole 
path of S. We impose the following assumptions, which will hold throughout: 

• Assumption (Al): All the coefficients are J^t-adapted and uniformly 
bounded in . 

• Assumption (A2): The volatility for the financial index P is uniformly 
elhptic: Hcrf (a;)|| > e > for aU (t.uj). 

• Assumption (A3): The payoff g is a positive bounded functional. 

The non-tradability of assets S could be applicable in many situations - it 
might be that these assets are (i) not traded at all, or (ii) that they are traded 
illiquidly, or (iii) that they are in fact liquidly traded but the investor concerned 
is not permitted to trade them for some reason. Our main application is to 
the counterparty risk of derivatives where the payoff depends upon both the 
value of the counterparty's assets and the asset underlying the derivative it- 
self. A second potential area of application is to residual or basis risks arising 
when the asset (s) used for hedging differ from the assets underlying the con- 
tract in question (see [13]). Typically this arises when the assets underlying 
the derivative are illiquidly traded (case (ii) above) and standardized futures 
contracts are used instead. Contracts may involve several assets, for example, 
a spread option with payoff (K — — S^)^ or a basket option with payoff 
(K — — ... — S^)^ . Such contracts frequently arise in applications to com- 
modity, energy, and weather derivatives. Finally, a one dimensional example of 
the situation in (iii) is that of employee stock options (see Henderson [H]). 

Our approach is to consider the utility indifference valuation for such con- 
tingent claims. For this we need to consider the optimization problem for the 
investor both with and without the option. The investor has initial wealth 
a; € M, and is able to trade the financial index with price Pt (and riskless bond 
with price 1). This will enable the investor to partially hedge the risks she is 
exposed to via her position in the claim. Depending on the context, the financial 
index may be a stock, commodity or currency index, for example. 

The holder of the option has an exponential utility function with respect to 
her terminal wealth: 

Urix) = for7>0. 

The investor holds A units of the claim, whose price is denoted as £g and 
is to be determine41, and invests her remaining wealth x — £q in the financial 
index P. The investor will follow an admissible trading strategy: 



G ^ad[0,T] = < TT : TT is /"(-adapted, and supi? 
I 

for any J't-stopping time r e [0,T].} 







J \7Tt\^dt 









< CO 



^For Cq we use the superscript A to emphasize the dependence of the price on the number 
of units held, and use the subscript to denote the price at time 0. If one unit is held, we 



simply write Co rather than Cj. 



5 



which results in the wealth: 



(2.3) 



The integrability condition for the trading strategies is essentially the BMO 
martingale property for nsdWs- However this condition is not restrictive if 
we only want to price and hedge contingent claims with bounded payoff, as the 
corresponding optimal trading strategy will satisfy this condition anyway. 

The investor will optimize over such strategies to choose an optimal n* by 
maximizing her expected terminal utility: 



sup E" 

Te^„<i[o,T] 



""(7r)+Ag(5.) 



(2.4) 



To define the utility indifference price for the option, we also need to consider 
the optimization problem for the investor without the option. This involves the 
investor investing only in the financial index itself. Her wealth equation is the 
same as (j2.3p but starts from initial wealth x and she will choose an optimal n* 
by maximizing: 



sup E'' 



(2.5) 



We note that (|2.5p is a special case of (|2.4p with A = 0. 

The utility indifference price for option is the cash amount that the investor 
is willing to pay such that she is no worse off in expected utility terms than 
she would have been without the option. For a general overview of utility 
indifference pricing, we refer to the recent monograph edited by Carmona [10] 
and especially the survey article by Henderson and Hobson [22- therein. 

Definition 2.1 (Utility indifference valuation of option and hedge) 

The utility indifference price o/A units of the derivative with payoff g{S.) 
is defined by the solution to 



sup E^ 

7veAo.dlO.T] 



7(X^ '-0(7r)+Ag(S.) 



sup E^ 



-ix^(^) 



(2.6) 



The hedging strategy for A units of the derivative is defined by the difference in 
the optimal trading strategies tt* — tt* . 

One of main features of the utility indifference price is nonlinearity, i.e. 
€q ^ ACq. Given the assumption of exponential utility, the price of the option 
and corresponding hedging strategy can be represented by the solution of a 
quadratic BSDE. 

Lemma 2.2 Suppose that Assumptions (Al) (A2) and (A3) are satisfied. Let 
(y, Z) be the unique solution of the quadratic BSDE: 



Yt^\g{S.) 



f{s,Zs)ds 



{Zs,dWs), 



(2.7) 
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where the driver f(s,Zs) is given by 




[A.g(5.)|J-t]. 



(2.8) 




The above type of quadratic BSDE (|2.7I) can be derived by standard Mar- 
tingale Optimality Principle (see, for example, Theorem 7 of Hu et al [53] and 
Section 3 of Ankirchner et al [1] in a Brownian motion setting, and Theorem 
13 of Mania and Schweizer [33] and Section 2.1 of Morlais [21] in a general 
semimartingale setting). In the Appendix, we provide a different proof for 
Lemma 12.21 where we do not use Martingale Optimality Principle. Instead we 
consider the problem from risk-sensitive control prospective, and employ Gir- 
sanov's transformation and comparison principle for quadratic BSDEs to derive 
(|2.7p . We should also mention that considering utility indifference pricing from 
risk-sensitive control prospective seems to be new. 

On the other hand, the well-posedness of quadratic BSDE (|2.7p is guaranteed 
by Kobylanski [21]. Indeed, by Assumption (A3) the terminal data Xg{-) is 
uniformly bounded. The driver f{t,z) for z = {z^, - ■ ■ , z'^) S K'^ is continuous 
in z, and moreover, by Assumptions (Al) and (A2) on the coefficients, f{t,z) 
satisfies 



Hence, by Kobylanski [3T], there exists a unique solution {Y,Z) to BSDE (j2.7p . 

The probabilistic representation formula (j2.8D means the utility indifference 
price can be written as a nonlinear expectation of payoff under the original phys- 
ical measure. However, this presentation is sometimes not convenient. First, the 
driver f{t, z) for z = (z^, • • • , z"^) G E'' is quadratic in z, which is in particular 
not Lipschitz continuous. Secondly, some nice properties for linear expectation 
such as Dominated Convergence Theorem may not hold anymore for nonlinear 
expectation. Two natural questions to ask are 

• Can the utility indifference price be written as a linear expectation of the 
payoff under some equivalent pricing measure? i.e. 



• Is the above probabilistic representation formula (|2.9p useful? 

In the remainder of the paper, we shall show both of these questions have 
positive answers. We first try to find the pricing measure Q in (j2.9p . Based on 
the solution {Y,Z) to BSDE (|2.7I) . we introduce the Doleans-Dade exponential 
in the following lemma: 



i/(i,z)i <c(i+iizin, 



and 



|V,/(t,z)| < C(l + ||z||), for t e [0,r], a.s.. 



(2.9) 
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Lemma 2.3 The Doleans-Dade exponential £{N) defined by 



£iN)=e^p^N~-[N,N]}, 



where 



Nt = - 



7 



iZs) 



2||^flP 

/or t G [0, T] is a uniformly integrable martingale. 

Proof. Wc first note that the solution Y of BSDE (|2.7|j is uniformly bounded 
and the process /g(Zs,(iyVs) is a P-BMO martingale whose proof can be 
found, for example, in Lemma 12 of Hu et al [23] . 

Next, we verify that the Doleans-Dade exponential £{N) is a uniformly 



: sup 



<Csup 



{Zt,dWt 



IZfW^dt- 



2||^flP 



Cdt\Tr 



2A^f 



e [0,T], 






2 







(2.10) 



Since /(^(Z3,dyVs) is a P-BMO martingale, i.e. 



sup£;^ 







2 








J {Zt,dWt) 




= sup£;^ 










r 





< oo 



for any J^f-stopping time r G [0,T], the first term in the RHS of (|2.10p is 
uniformly bounded. Clearly, the second term is bounded by C^T. Hence N 
is a P-BMO martingale, and the Doleans-Dade exponential £{N) is uniformly 
integrable. ■ 

Since the Doleans-Dade exponential £{N) is uniformly integrable, we can 
change the probability measure from P to Q by |g = £:(7V). Under the measure 
Q, BSDE (1121) reduces to 



Yt = Xg{S.)-h I f{s,Zs)ds 

i-T 



(Zs,dWs 



Xg{S.) 



{Z„dBs 



where B = {B^, ■ ■ ■ , B'^) defined by S = W — [W, N] is the Brownian motion 
under Q by Girsanov's theorem. In other words, we have derived 

rt = i?Q[A5(5.)|J-t]. 



^We recall a continuous martingale M with [M, Af] < oo is called a P-BMO martingale 

if 

supE^[\MT - MtI^IJ"^] < oo 

T 

for any J-t-stopping time t S [0,T]. If M is a P-BMO martingale, its Dolc5ans-Dadc exponen- 
tial £{M) is in Doob's class "D, and therefore uniformly integrable. 
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Theorem 2.4 (Pseudo linear pricing rule for utility indifference valuation) 

Suppose that Assumptions (Al) (A2) and (A3) are satisfied. Then there 
exists a pricing measure Q defined by ^ = £{N), where the Doleans-Dade 
exponential £{N) is defined in Lemma \2.3l such that the utility indifference 
price Cq can be written as the linear expectation of the payoff g under Q ; 

= E^[Xg{S.)]. 

Proof. The main steps are already shown above. What is left is that 
we only need to show the martingale representation part Z is invariant under 
the equivalent change of probability measure. Actually, this follows by the 
uniqueness of the special semimartingale decomposition, whose proof can be 
found in Lemma 7 of Liang et al [37^ . ■ 

The probabilistic representation of the utility indifference price (j2.9p appears 
to be "linear" at first glance. Notwithstanding, the nonlinearity is hidden in 
the pricing measure Q, which depends on the price £q through the Doleans- 
Dade exponential £{N). Indeed, if the unit indifference price is given by Cq = 

[g{S.)], where Q is defined by the Doleans-Dade exponential which depends 
on the price Cq, then, 

€^ = i?Q[A.9(5.)] = XE%iS.)] ^ XE^' [g{S.)] ^ XU- 

We are not the first to discover such a pseudo linear pricing rule for the 
exponential utility indifference price. A similar probabilistic representation to 
Theorem [23] is given in Proposition 11 of Mania and Schweizer [W, where their 
payoff is a general random variable. In contrast, since we specify the dynamics 
of the underlying assets and the payoff structure, a fully coupled linear FBSDE 
appears naturally. Indeed, if we write down the equations of the non-traded 
assets S under the measure Q: 

f ^{'•:-|<"^'^'> -lit -f)}^- 

together with the dynamic of the utility indifference price under the measure 

Yt^Xg{S:)^J^ {Zs,dBs), 

we derive a fully coupled linear FBSDE. Since the associated equations are linear 
(but with a loop) , Picard iteration can be employed to approximate the solution 
and the corresponding utility indifference price, which results in the following 
approximation, whose proof is postponed to the Appendix. 

Theorem 2.5 (Linear approximation for utility indifference price) 

Suppose that Assumptions (Al) (A2) and (A3) are satisfied, and moreover, 
for either T or K small enough, the payoff functional satisfies 

n 

\g{S.) - g{S.)\ <K sup ^ | \nSl - \nSl\. 
te[o,T] 
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Let B = {B^, • • • , B"^) be the Brownian motion under a given probability mea- 
sure Q. Define the fo 
= (s, • • • , s), and 



sure Q. Define the following sequence {5™ = (S"™'^, • • • , S'^ ")}m>a iteratively: 



^ Mt - 7;\^t^^t I + ON PII2 v'^t ' A ) ]>dt+ {at,dBt), 



2' *' * ' 2||CTf||2 

where Z,™ = (Z™'^, ■ ■ ■ , Z™'"^) is the martingale representation of Xg{S, 

Then the utility indifference price Cq is approximated by 

= E^iXgiSni 

and 

lim c^'^^g:^. 

m— >-oo 



To summarize, Theorem 12 .41 and Lemma l2 . 2 1 provide two alternative ways to 
represent the utihty indifference price: 

€^^E^[Xg{S.)]=£^[Xg{S.)]. 



One apphcation of our representation in Theorem l2 .41 was illustrated in Theorem 
12.51 A second application will be given in Section [3] where we shall utilize the 
representation to prove the convergence of semigroup approximation for the 
utility indifference price by splitting method in Theorem 13.81 



3 The Markovian Case and the SpUtting Method 

In this section, we specialize to a Markovian setting by assuming that all 
the coefficients are constant, denoted by fip, ap, fj-i, at etc., and the payoff g is 
a function rather than a functional. 

• Assumption (Al)': All the coefficients are nonzero constants. 

• Assumption (A3)': The payoff g is a positive bounded function. 

We consider a one factor model where the prices of the traded financial 
index P and each of the non-traded assets S** are driven by one common market 
Brownian motion, as well as an independent Brownian motion, representing the 
idiosyncratic risk of each asset. In other words, d — n + 2, erf"' — ii j n + 1 
and n + 2, and a^' = if i 7^ j and j ^ n + 1. The price processes become 

dP 

= Updt + apdWl'+^ + CTp^iW^i"+^ (3.1) 

" t 

and 

—f = fidt + adWl + adW^+\ (3.2) 
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To simplify the notation, we set ap = crp^n+ii = (7p,n+2, = o'u, and 

The price of each non-traded asset reflects exposure to the traded or 
market risk Ty+i through volatihty ai and idiosyncratic risk through id- 
iosyncratic volatility ctj. We define the following parameters for the financial 
index P: 



2 ' 



„P _ A^PgP . aP _ 

{opr 





{dp 


2 + (ap)2 




^pdp 


{ap 


? + i^P? 






{ap 


Y + {^p? 



{^pY 

We first characterize the utility indifference price by the solution of the 
following semilinear PDE with quadratic gradients. 

Lemma 3.1 Suppose that Assumptions (Al)' (A2) and (A3)' are satisfied, and 
moreover, the payoff function g(-) is Lipschitz continuous. Then the utility 
indifference price of X units of the option with payoff g {St) is given by £^(-,0), 
where £^ is the classical solution of the PDE: 



-i E", = l [(1 - S^SjdsX^ds^<t^ = 0, 

^\;T)^Xg{-) 



(3.3) 



on the domain {s,t) G M" x [0,T], where s — (si,--- ,s„), and the hedging 
strategy for A units of the option is given by 

_p n 

-'^y2a,s,dsX\ (3.4) 

Since the proof of Lemma 13.11 follows from the standard Dynamic Program- 
ming Principle, we omit it (for example see |18j in one dimensional setting) . We 
only remark that the classical solution to p.ip exists, and moreover, and its 
gradients are uniformly bounded (see, for example. Theorem 2.9 of Delarue [T^ 
for the proof). We also note that the number of units A only appears in the 
terminal condition. In the following, we present the case A = 1, and the price 
is simply denoted by £. 

In general, there is no explicit solution to PDE p.3p . However, if the traded 
financial index is independent of the non-traded assets, i.e. ap = 0, the explicit 
solution is available by the Cole-Hopf transformation (see comments on the one- 
dimensional case after Proposition 13. 6p . Essentially, this means the financial 
index P is not useful as a hedging tool for the investor. 
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Corollary 3.2 Suppose that ap = 0, and Assumptions (Al)' (A2) and (A3)' 
are satisfied, and moreover, g{-) is Lipschitz continuous. Then the utility indif- 
ference price of the option with payoff g [St) is given by 

1 



7 



g-7S('ST) 



(3.5) 



We now contrast the above to the situation if the market were complete. 
If the underlying assets S = (S*^, • • • , S'") could be traded, the market would 
become complete, and the pricing and hedging of the contingent claim with 
payoff g(5T) falls into the classical multidimensional Black-Scholes framework. 

Corollary 3.3 Suppose that 5 = (S'\ • • • ,5'") are traded assets, and Assump- 
tions (Al)' (A2) and (AS)' are satisfied, and moreover, g{-) is Lipschitz con- 
tinuous. Then the price of the option with payoff g [St) is given by £, where £ 
is the classical solution of the PDE: 

dt€{s,t) + ^J2'Li ^1^^190,8,^+ ^J2to=i^^^jS,Sjds,sj€ = 0, 

(3.6) 

£(-,r)=g(.) 

on the domain {s,t) G R" x [0,T]. 

Based on the pricing equation p.3p . we will present a number of properties 
of the option price. First notice that the pricing equation p. 31) has additional 
nonlinear terms of quadratic gradients relative to the complete market PDE in 

am . 

Our first property of the utility indifference price concerns risk aversion. 
Intuitively, the more risk averse the option holder is, the less she would be 
willing to pay for the option. 

Proposition 3.4 // the risk aversion parameter 7 increases, the unit utility 
indifference price € will decrease. 

Proof. The proof is based on the comparison principle for PDE p. 31) . 
Suppose < 71 < 72. The corresponding PDEs are denoted by PDE'^^ and 
PDE^^ respectively, and their solutions are denoted by and respectively. 
Since the terms involving 7 can be regrouped as 

n n 
i=l i=\ 

SO all the coefficients of 7 are less than zero, we have 

PDE'^^{V-) > 0. 

Hence is the subsolution to PDE'^^ . On the other hand (t^'^ is the super- 
solution to PDE'^^, and €''^\t=T — 't^'^\t=T- By the comparison principle, we 
conclude < on x [0,T]. ■ 

Proposition 3.5 Assume that 

= (3.7) 

for i = 1, • • • ,n. Then the unit utility indifference price £ is decreasing in the 
idiosyncratic volatility of the traded asset <7p (or its proportion of total volatility, 
1-rP). 
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Proof. The proof is again based on the comparison principle for p.3p . and 
is similar to that of Proposition !?^ ■ 

This tells us that the higher the idiosyncratic volatility dp of the traded asset 
(or as a proportion of total volatility) , the worse it is as a hedging instrument, 
and the lower the price one is willing to pay. Proposition 13.51 generalizes the 
monotonicity obtained in the one dimensional non-traded asset model (see for 
example, Henderson |19) and Frei and Schweizer |16) in a non-Markovian model 
with stochastic correlation). 

The restriction p.7p in fact corresponds to a relation between the Sharpe 
ratios of the non-traded assets S and the financial index P. Define the Sharpe 
ratio of 5" to be SRi = Hi/ \/ of + g1 and similarly, the Sharpe ratio for the 
financial index P by SRp — [ipj \Jo'\y + a\. Then (13. 7p is equivalent to the 
relation 

SRi — I — , = — — I SRp = PipSRp 

where pip is the correlation between 5' and P. This corresponds to the relation 
we expect from the capital asset pricing model (CAPM) when assets are traded. 
Since not all assets are traded here, we would not necessarily expect p.7p to 
hold. 

The final result in this section concerns recovery of the complete market 
price given in Corollarv 13.31 

Proposition 3.6 (Asymptotic results) 

• (i) If — ^ for i — 1, ■ ■ • ,n, then the unit utility indifference price £ 
will converge to the complete market price £ as 7 — >■ 0. 

• (ii) If ^ — ^ for i = I,-- - ,n, then as ap,ai — > the unit utility 
indifference price £ will converge to the complete market price £ with 

CTp, fJi = 0. 

Proof. We prove (i), whilst the proof of (ii) is similar. Under the assump- 
tions in (i), when 7 — > 0, by the Arzela-Ascoli compactness criterion, there 
exists a subsequence 7„ — > such that the solutions of (j3.3p . denoted by C"", 
uniformly converge to £ on any compact subset of M" x [0, T], where £ satisfies 

daH). ■ 

Again, we see in Proposition 13.61 that the CAPM restrictions must apply to 
the Sharpe ratios. The intuition here is that when the idiosyncratic volatilities 
disappear, and when assets are traded, there cannot be a difference in using the 
financial index P or the assets themselves to hedge. 

For a semilinear PDF with quadratic gradients like (j3.3p . it is not usually 
possible to obtain an explicit solution. A special case where an explicit solution 
does exist is the one dimensional version. Taking n = 1, d — 2 and api — 
in p.3p recovers the pricing PDF of [H], [TH] and [33], which is solved by the 
Cole-Hopf transformation. However, this transformation does not apply directly 
to our multidimensional problem (13.31) because the coefhcients of the quadratic 
gradient terms in (j3.3p do not match. We note that even applying standard finite 
difference method to numerically solve the PDF p.3p is troublesome due to the 
high dimension and nonlinearity. Instead, we will develop a splitting algorithm 
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which will enable us to take advantage again of the Cole-Hopf transformation 
to linearize the PDEs. 

Splitting methods (or fractional step, prediction and correction) can be dated 
back to Marchuk [?D] in the late 1960's (see also [2]). Application of splitting 
idea to nonlinear PDEs such as HJB equations is difficult mainly because of the 
verification of the convergence for the approximation scheme. This difficulty was 
overcome by Barles and Souganidis [1], who employed the idea of the viscosity 
solution and proved that any monotone, stable and consistent numerical scheme 
converges (to the correct solution) provided there exists a comparison principle 
for the limiting equation. Barles and Souganidis [4] pointed out that their result 
could be used to justify most standard splitting methods. 

The idea of splitting in our setting is the following. The local time derivative 
of the pricing PDE p. 31) depends on the sum of semigroup operators (or the 
associated infinitesimal operators) corresponding to the different factors. These 
semigroups usually are of different nature. For each sub-problem corresponding 
to each semigroup there might be an effective way providing solutions. For the 
sum of these semigroups, however, we usually can not find an accurate method. 
Hence, application of splitting method means that instead of the sum, we treat 
the semigroup operators separately. 

Of course, the tricky part is how to split the equation (or how to group 
factors) effectively. In the following, we separate the pricing PDE (13. 3|) into two 
pricing factors. Take A = 1. We first make the log-transformation: Xi = Insi, 
and define a new operator: 

dri ~* dxi ^ 

i—l 

Then p.3p reduces to 

-. n n n 

%—\ i—\ i—\ 

+ i5^r,e:- J(i-s^)(a^e:)2 = o, (3.9) 



where 



A,=y.,-h^al+a\)-¥a,. (3.10) 



Define two operators: 



2—1 i=l i—l 

For any < Ti < Ta < T, and any smooth function (j) e C°°(IR™) for m = 
1, n, we define the following nonlinear backward semigroup operators S'(T'i, T2) 
C°^(K™) ^C°°(R™) by 0(-) h-> e(-,ri) where 

dt€' + h'e = o; £■'(•, T2) = (/)(•) 

on the domain [Ti, T2] x R™ for i = 1, 2. 



14 



We observe that one of the equations can be hnearized by a Cole-Hopf trans- 
formation. Indeed, by letting £^ = exp(— 7(1 — k^)€^), then we have £^ satis- 
fying 

9tl' + ^dnnl' = 0. (3.11) 
By letting £^ = exp(— 7^^), then we have satisfying 

i=l i=l 

Lemma 3.7 The operators S'(Ti, T2) /or i = 1,2 have the following properties: 

• (i) For any smooth function ^ G C°°(]R'"), 

\imSHTuT2)ct> = 4> 

-I lT-'2 

uniformly on any compact subset o/R"*. 

• (a) For any < Ti < T2 KT^ < T, 

S\T,,T3)cf> = S\T,,T2)S\T2,Ts)cf>. 

• (Hi) 

S\T2,T2)ct> = c/,. 

((i) (a) and (Hi) ensure that S*(Ti,T2) is indeed a semigroup operator.) 

• (iv) If (j) > 'ij) where tp is another smooth function, then 

s'(ri,r2)<^> s'(Ti,r2)v. 

• (v) For any constant k gM., 

s\TuT2m-)+k) = s\n,T2m-))+k. 

• (vi) 



TitT2 T2-T1 

uniformly on any compact subset o/M"*. 

Proof, (i)-(v) are immediate. We only prove (vi) in the following. We first 
prove the case i = 1. Note that S^{Ti,T2) = 8^(0, T2 - Ti), and it is enough to 
prove that for any (f) e C°°(R), 

lim 

T2-T^i0 T2-T1 

uniformly on any compact subset of R. 
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Let {^l,T,Tt,P) be a filtered probability space satisfying the usual condi- 
tions, on which supports Brownian motion W with P(Wo = tj) = 1. By Ito's 
formula, we have 



T2-Ti 



T2-T1) - 



T2-Ti 



-{l-R''){d^€'f{s,W,)ds 



T2-T1 



d^€\s,Ws}dWs 



for t e [0, T2 — Ti] . We note that dr,€^ {t, rj) is uniformly bounded, and therefore, 

Nt = - f J(l - K^)a^£i(s, W,)dWs, for t e [0,r2 - Ti] 
Jo ^ 

is a P-BMO martingale, which implies that the Doleans-Dade exponential £{N) 
is uniformly integrable. Hence, we can define a new probability measure Q by 
^ = £(N). By Girsanov's theorem, 

Bt = Wt- [W, N]t = Wt+ f J (1 - K'')dr,C^ (s, W,)ds, for t e [0, T] 

Jo ^ 

is Brownian motion under the probability measure Q, and moreover, 
€\t,Wt) = (f>{WT2-T,) - p dr,<L\s,Ws)dB, 

= E'^[<t>{WT2-Tm. 

Therefore, by Ito's formula, we have 

s\Q,T2-T^)m-m 
T2-T1 

' dr,<t>{Ws)dWs + ' ^dr,r,^{Ws)d{W) 
T2-T1 /-Ta-Ti 



T2 






1 


T2 


-Ti 




1 


T2 






pT2 — 


+ 


Jo 




1 


T2 





dncP{Ws)dB, - [ ' ' J(l - K'')dr,(P{Ws)dr,<i\s,W,)ds 
Jo ^ 



i 

2 



dr,n(l){Ws)ds 



T2-T1 



Note that 9^£^(0,ry) = dn(j){r}) and by Dominated Convergence Theorem, 
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The proof for the case i = 2 is similar, so we only sketch its proof. We apply 
Ito's formula to €'^{t, Xt), where X = (X^, ■ ■ ■ , X^) is given by 

Xi = x' + [ Aids + / adWl 
Jq Jq 

Therefore, by changing the probability measure, we obtain 

^=l 

= E'^[cj,{XT,^T,)\:Ft], 
where Q is defined by the the Doleans-Dade exponential £{N) with 

Nt = - l'^^9x,^^is: Xs)dWl, for t G [0,T2 - T,]. 

•^0 1=1 

The rest of the proof follows by the same argument as for the case i = 1. ■ 

Next we use semigroup operators S^(Ti, T2) and S^(Ti, T2) to give the semi- 
group approximation for the solution of PDE p.Op for PDE p.3p ). which is the 
main result of this section. 

Theorem 3.8 (Semigroup approximation for utility indifference price) 

Suppose that Assumptions (Al)' (A2) and (A3)' are satisfied, and moreover, 
the payoff function g{-) is Lipschitz continuous. Let tt : = to < ti < ■ ■ ■ < 
tN = T be the partition of [0, T] with mesh: 

|7r| := max \ti+i - ti\. 

0<i<N~l 

Then the unit utility indifference price of the derivative with payoff gi^Sr), which 
is denoted by £(-,0) is approximated by 

e{;0)^ n S\t,,t,+,)S'(t,,t,+i)g{-), 

and 

lim £^(-,0) = £(-,0). 
uniformly on any compact subset o/R". 

Proof. The proof is based on the Barles-Souganidis monotone scheme @], in 
which they proved that any monotone, stable and consistent numerical scheme 
converges (to the correct solution) provided there exists a comparison principle 
for the limiting equation. 

In the following, we verify the above conditions. For any < Ti < T2 < T 
and any smooth functions (/> > "0? by (iv) in Lemma 13.71 

Si(ri,T2)S2(Ti,r2)0 > S\Ti,T2)S\Ti,T2yj, 
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so the scheme is stable. 

Finally, we verify the scheme is consistent: 





+ Lis2(Ti,r2)0 




-Lis2(ri,r2)<^ + LV 
= (/) + (//) + (///). 



(3.13) 



By (vi) in Lemma 13.71 the terms (7) and (//) in (j3.13l) converge to when 
Ti t 72- By (i) in Lemma [?771 the term (///) in (|3.13p also converges to when 
Ti t 22. Therefore the scheme is consistent. ■ 

4 Application to Counterparty Risk of Deriva- 
tives 

In this section, we apply our multidimensional non-traded assets model (in 
particular, the one factor model of Section 3) to consider the counterparty risk 
of derivatives. Our concern as the buyer or holder of the option is that the 
writer or counterparty may default on the option with payoff h{S.) at maturity 
T and we will not receive the full payoff. We have in mind several examples. A 
natural example is that of a commodity producer who is writing options as part 
of a hedging program (eg. collars). Some of these options may be written on 
illiquidly traded assets and thus the option holder is subject to basis risk and 
in addition, is concerned with the default risk of the option writer. A second 
example is the default risk of a financial institution who has sold options on 
various underlying assets - stocks, foreign exchange or commodities. In addition 
to the possibility of basis risk, the buyer of these options does not always have 
the ability to trade the underlying asset, or perhaps they choose not to (they 
may be using the derivative as part of a hedge already). A further example 
may be that of a purchaser of insurance concerned with the default risk of the 
insurer. Typically the insured party does not trade at all, which motivates 
our consideration of this special case. Finally, the option holder may be an 
employee of a company who receives employee stock options if the company 
remains solvent. She is restricted from trading the stock of the company, but can 
trade other indices or stocks in the market. In contrast to the other examples, 
here the assets of the counterparty and the underlying stock are those of the 
same company. 

We consider an option written on an underlying asset with price with 
payoff h{S)p) at maturity T. Counterparty default is modeled by comparing the 
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value of the counterparty's assets S*^ to a default threshold D at maturity, which 
depends on the liabilities of the counterparty. Following Klein we consider 
the situation D = L, where L refers to the option writer's liabilities, assumed 
to be a constant]! If the writer defaults, the holder will receive the proportion 
h{S^)/L of the assets S"^ that her option represents of the writer's liabilities, 
scaled to reflect a proportional deadweight loss of a e [0, 1]. The payoff of the 
vulnerable option taking counterparty default into account is 

giS'r,S'r) = h{S'r)lis-^>L} + (1 - a)i^S'rhs^^<L}- (4.1) 

To guarantee g is bounded, a sufficient condition is that the payoff h{S^) is 
bounded, for example, a put option. Note that there is a singular point of g at 

— L, so we have to approximate g by a sequence of (nondecreasing) Lipschitz 
continuous functions g'^. For the numerical simulation, we only need to choose 
one g"^ for e small enough. 

The underlying asset and the value of the counterparty's assets are 
both taken to be non-traded assets so ti = 2 and prices follow (13. 2p . The option 
holder faces some unhedgeable price risk (due to S^) and some unhedgeable 
counterparty default risk (due to S"^). She can partially hedge risks by trading 
the financial index P following (13. ip . 

Our benchmark models for comparison are those of Johnson and Stulz [2 8) , 
Klein [53] and Klein and Inglis [5^ who take a structural approach to price 
vulnerable options in a complete market setting and obtain two dimensional 
Black Scholes style formulas. Other treatments of vulnerable options include 
complete market models with intertemporal default (see Liang and Ren |38| and 
the references therein) . Implicit in this prior literature are the twin assumptions 
that the asset underlying the option and the assets of the counterparty can be 
traded, and therefore, can be used to hedge the counterparty risk of derivatives. 
Our use of the utility indifference approach is motivated by its recent use in 
credit risk modeling where the concern is the default of the reference name 
rather than the default of the counterpartylf] Several authors have applied it 
in modeling of defaultable bonds where the problem remains one dimensional, 
see in particular Leung et al [34] and Liang and Jiang [35] and Jaimungal and 
Sigloch [25]. In contrast, options subject to counterparty risk are a natural 
situation where two or more dimensions arise. 

We were also motivated to study indifference pricing of derivatives subject 
to counterparty risk by the work of Hung and Liu [23] and Murgoci [33] • These 
papers take a good deal bounds approach to pricing, acknowledging the incom- 
pleteness of the market but producing prices which are linear in quantity. Fur- 
thermore, the method does not allow for any partial hedging on the part of the 
investor and can produce bounds which can be quite wide. 

Based on Theorem 13.81 we give the following approximation scheme for the 
unit utility indifference price £(si,S2,0) of the vulnerable option. Following 
p.8p we make the change of variable x = Insi and y = lns2, and define a new 

^Generalizations to D = f{S^) are easily incorporated and allow for the option liability 
itself to influence default, eg. f{x) = h{x) + L was considered by Klein and Inglis '30' in a 
risk neutral setting. 

^ Utility based pricing has also been utilized by Bielecki and Jeanblanc 7,, Sircar and 
Zariphopolou 48 and recently Jiao et al 1261 |27l in an intensity based setting. 
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operator: 

d _ d _ d 
or] ox ay 

• (i) Partition [0,T] into N equal intervals: 

Q = ta<ti<--- <tN =T. 

• (ii) On [iAr_i,tjv], predict the solution by solving the following PDE with 
the given terminal data g'^: 

dt^ + i^^^ci - ^(1 - s^)(9^g:1)2 = 0, 

The above equation can be linearized via the Cole-Hopf transformation: 

£i = exp(-7(l - fi^)€^). 
In particular, we obtain (i^\t=tN-i- 

• (iii) On [tN-i,tN], correct the solution by solving the following PDE with 
the terminal data (t^\t=tN-i- 

where Ai,A2 are given in p.lOp to be Ai — fii — ^{af + (jf) — 'd^(Ji; i = 1,2. 
The above equation can also be linearized by making the exponential 
transformation: 

= exp(-7e:^). 

In particular, we obtain £^|t=t„_j, which is used as the approximation of 

• (iv) Repeat the above procedure on ^Af-i], and obtain £|t=t„_2 ' ' ' 

We present results for the European put with payoff h{S^) = {K — S^)^ . If 
and S*^ are positively correlated, this means when the put option is valuable 
(in-the-money), the firm's assets S*^ tend to be small, so there is a high risk 
of default. It is important to take counterparty risk into account for puts in 
this case, as it will have a relatively large impact on the price. (This would 
be even more significant when the default trigger involves the option liability). 
However, for a call, when the call is in-the-money, there is little default risk, and 
so counterparty risk is less important. Unless otherwise stated, the parameters 
are: K = 150; T = I; = 50; 5^ = 100; L = 1000; a = 0.05; = 1: fip = 0.1; 
crp = 0.15; ap = 0.2; = 0.15; ai = 0.25; cti = 0.3; /za = 0.1; CT2 = 0.3; 
(72 — 0.2. These parameters result in correlation between the underlying asset 
and firm's assets of pi2 = 0.4; and correlations between each asset and the 
financial index P of pip = 0.6; p2P = 0.4. . 
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In Figure [T] we show how the approximation converges as we increase the 
number of time steps N. For our parameter values, A'^ = 11 steps is sufficient 
for the prices to converge and we use it in all subsequent figures. We aim to 
compare the utility indifference price with hedging in the financial index with 
the benchmark risk neutral price in a complete market (computed as in Lemma 
13.31 with n — 2, as studied in Klein [12 )• We also compare to the situation 
where the financial index is independent of the other assets and thus there is 
no hedging carried out. This price was given earlier in Corollary 13.21 Figure [5] 
provides a demonstration of the accuracy of the algorithm. We take tip = 
and compare the splitting approximation to the formula in Corollary 13.21 

Figure [3] shows the vulnerable option price(s) against the underlying asset 
price S^. The two panels of Figure [3] are intended to illustrate a "close or likely 
to default" scenario (the left panel with = 500 relative to L = 1000) and a 
"far or unlikely to default" scenario (the right panel with S2 ~ 1400 relative to 
L — 1000). As expected, in both panels the risk neutral or complete market 
price is the highest. In both panels, as the underlying asset price becomes 
very large, all option prices tend to zero, as the put is worthless, regardless of 
the default. At = 0, in the right panel, the option price is equal to the 
option strike K = 150. In the left panel, the option price is lower due to the 
risk of counterparty default. As increases, all option prices decrease, as the 
moneyness of the put changes. When is close to zero, we see a dramatic 
drop in the utility indifference prices (relative to the risk neutral prices) due 
to the risk aversion towards unhedgeable price and default risks. Recall that 
since the underlying asset and firm's assets are positively correlated, default 
risk becomes more important for low values of the underlying asset. The price 
drop is much more significant in the left panel (and in this extreme case, the 
option price drops down to zero if no hedging can be carried out), where the 
likelihood of counterparty default is higher. The option holder's risk aversion 
causes the utility indifference prices to lie below the risk neutral price (in each 
default scenario) with the relative discount to the risk neutral price being much 
greater in the left panel where default is more likely. Assuming the holder can 
hedge in the financial index, there is a drop of around 75% from the risk-neutral 
price to the utility indifference price. In the right panel, where the likelihood of 
default is relatively low, the difference between the utility indifference price(s) 
and the risk-neutral price is not as dramatic, and is at most around 20% of 
the risk-neutral price. We also see that the ability to hedge in the correlated 
financial index (versus no hedging at all) is more important when the default 
risk is higher (in the left panel). 

Figure [3] displays the impact of the option writer's asset value on the 
option price for a fixed asset price S^. We see a dramatic difference in the 
behavior of the risk neutral price and the utility indifference prices. Under risk 
neutrality, the option price increases smoothly with S"^. However, under utility 
indifference, the prices are low and do not change much with values of 5^ below 
the default trigger of L = 1000. This is despite the put being in-the-money. As 

increases beyond the default level, the likelihood of default diminishes, and 
the utility indifference prices start to increase with S'^. Note that the utility 
price is not always below the risk neutral price. Although ProDOsition l3.6l tells us 
that the risk neutral price is obtained as a limiting case of the utility indifference 
price, it requires condition p.7p to hold. 

Figure compares how the vulnerable option price changes with risk aver- 
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sion parameter 7, and the idiosyncratic volatilities cri,cr2-Tlie left panel plots 
vulnerable option prices against maturity T for various values of risk aversion 
7. We see that the more risk averse the option holder is, the less she will pay 
for the option, consistent with Proposition 13.41 The other observation is that 
option prices for a fixed 7 are decreasing with maturity T . The risk neutral 
price is also decreasing with T, albeit very gradually. This is in contrast to risk 
neutral prices for non-default European put options which will increase in T 
provided there are no dividends. The reason is that there is a tradeoff between 
price and default risk. If the maturity is longer, there is more chance for both 

and to fall - falling means the put is more valuable, but S'^ falling 
increases the default risk. For the parameters considered, the default risk is the 
dominant factor and thus the option price decreases with T . This is also in 
contrast to the call option, where Klein [55] reports that the risk neutral price 
increases with maturity. 

Recall that we do not expect price monotonicity in terms of the correlations, 
except in the situation outlined in Proposition [3?5] Here we give an example of 
prices for various values of the idiosyncratic volatilities ci , (T2 ■ The left panel 
sets parameters to be ni = 0.1 and ^.2 = 0.06 to satisfy the CAPM restriction on 
Sharpe ratios. If cti = 0-2 = 0, then we have pi2 = 1, pip = 0.8 and p2P — 0.8. 
Similarly if cti — 0.25 and (T2 = 0.3 then pi2 = 0.4, pip = 0.6 and p2P = 0.4. 
Finally, if cti = (T2 = 1 then pi2 = 0.06, pip ~ 0.2 and p2p = 0.16. We see 
that as (Ti,cr2 increase, the utility indifference price falls. Correspondingly, as 
the correlations P12, Pip, P2P increase, the option price rises. 

5 Concluding Remarks 

We close with several suggestions for further work. In this paper, we prove 
the convergence of our splitting method for the utility indifference pricing PDE 
but do not concern ourselves with the convergence rate. Tools such as Krylov's 
idea of "shaking the coefficients" (see [23]) and Barles and Jakobsen's "switching 
system" (see [21) could be useful. A second suggestion for further investigation 
is to employ the splitting method to give a numerical scheme to solve quadratic 
BSDEs. To achieve this goal, the first problem may be to understand the mono- 
tone scheme from a probabilistic perspective, we refer to Cont and Voltchkova 
[g and Dai and Wu [12 . 

Several developments in the Markovian setting of Section |2| may be of inter- 
est. A first comment is that the splitting relies on the special setup of the one 
factor model and it does not seem possible to extend beyond this case. It would 
be interesting to consider alternative numerical methods which could be used in 
the general multidimensional framework. Second, following Eberlein and Madan 
[15] we could extend our counterparty risk model in Section U) to treat assets 
and liabilities as separate stochastic processes. Finally, in this paper we price 
European claims in a multidimensional non-traded assets framework. Exten- 
sions to corresponding American claims in a multidimensional model (following 
work of Oberman and Zariphopoulou [46| and Henderson [20) in one dimension) 
would allow us to consider intertemporal counterparty default. 
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Appendix 



Proof of Lemma 12.21 We characterize the vahie function of the optimiza- 
tion problem (12.41) by the solution of quadratic BSDE. By plugging (|2.3p into 
(|2.4p . we have 



sup E 

Te^„d[o,T] 

— e 



.-7 (/o^ 'r, (^f (af .rfW, ) ) + Ag(5. ) ) 



exp<-7 sup 3^0 (tt) 



where 3^o(i') denotes the risk-sensitive control criterion: 



-1 



3^0 W = —lnE^ 

1 



We further introduce the risk-sensitive control problem: 



3^0 = sup 3^0 (tt)- 

In the following, we characterize both yo{T:) and 3^o by the solutions of quadratic 
BSDEs. 

Let us first consider the risk-sensitive control criterion 3^o (■"') ■ For any given 
trading strategy tt € Aad[0,T], we define P-BMO martingale: 



77r,(af,dW,), forie[0,r]. 



Indeed, for any /"(-stopping time r € [0,T], we have 



J%'\\afr\nA'ds 



supi;^[|A^T - A^rH-T^r] = sup 

r T 

< CT sup 



Hence the Doleans-Dade exponential is uniformly integrable, and we 

change the probability measure by defining %^£{N). Under the new proba- 
bility measure Q, the risk-sensitive control criterion becomes 













< OO 









3^0 W = — InSQ 

7 



which is the unique solution to the following quadratic BSDE: 

yt{TT) ^ Xg{s.) + (^^,W-|||z,||2)ds~ / {Zs,dm) (0.1) 



with 



2|I~-P||2|_ |2 
2' 
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and VV = W — [W, AT] being the Brownian motion under Q. Indeed, note that 
(jO.ip can be reformulated as 

t T T T 

yt{^)+ I Fs{TT)ds ^ Xg{S.) + [ Fs{n)ds- [ ^\\Zs\\^ds- [ 

Jo JQ Jt ^ Jt 

By changing variables: 

we have ^ 
which has the explicit solution: 



,-7(Ag(5-)+/o^i=^.Wds) 



e-/o^(wf^=-^7'lkf||^|7r,|^)dSg-7Ag(5.)' 



Next, we consider the risk-sensitive control problem for 3^o- Note that under 
the probability measure P, (|0.ip becomes 



l\\Zs\? 
2" " 



{Zs^dWs 



(0.2) 

By the comparison principle for quadratic BSDE, 3^o — sup^g_^^^[o t] 3^o(7'') is 
the unique solution to the following quadratic BSDE: 



yt^\g{S.) + J^ Fsds-J^ (Z, 



dWs 



with 



=sup{^^,W -7«,Z,)7r, - 



2" '^'l 2||af||2 



7 



and 



To verify that tt* € ^Qd[0, T], we only need to note that Jg{Zs,d'Ws) is a P-BMO 
martingale and Assumptions (Al) and (A2) on the coefficients. 

The optimization problem (j2.5p is a special case of (12.41) with A = 0, whose 
solution is 

/ ( 1 

— exp < —7 X 



27 Jo Ik 



P||2 



and the optimal trading strategy tt* is given by tt^* 
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Finally, by Definition 12. 1[ the price €q is given by the solution to 



)+As(5.) 



2t Jo ji„-P||2 



Therefore, 



T |,,P|2 



27^0 \W 



P\\2 



ds 



Jo Jo ^1 Jo iWiW 



ds 



Xg{S.)+ / /(s,Z,)ds 



Zs-dWs, 



and the hedging strategy for such A units of the option is given by 



P||2 



which completes the proof. 



Proof of Theorem 12.51 The type of FBSDE is in fact a special case of 
functional differential equations studied by Liang et al in [36^ and [37 . We define 
©([0,T],R"), the space of continuous and J^t-adapted processes valued in R" 
such that supjg|o_7-] J27=i \'^t \ G L'^{^t^t,Q) and endowed with the norm: 



1/2 



||5||6=£;^ sup J2\si\' 



Then for 5™ and S"^'^^, we have 

E^^l sup y lln^r^^'^-ln^r'f 1 

it6[0,T]-^ J 

= E^\ sup ± fl{al, 
[te[o,T] Jo ^ 

< cVte"^ i^J^ \\z^ 



rrn rim 



1/2 



-1\ 7(0-^ erf) 



P 'T-m 'T-m— 1 



)ds 



1\ 1/2 



cVte^ 



\dBs) 



2-11/2 



1/2 



<CXVr£;'*| sup y llnS^'-ln^r"' 



1/2 



i|2 
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We iterate the above inequality and obtain 

||ln5"'+i-ln5™||e < (C/^\/T)"|| In^^ - ln5°||e. 
Hence if either T or if is smah enough, for any natural number p, 

p 

llln^^+f -ln5™||e < ^ || In5'"+J' - In5"+J'~i||e 



~ 1-ckVt 

when m oo. Therefore, In^™ is a Cauchy sequence in (3([0,T];R") and 
converges to some In 5, and moreover, by Dominated Convergence Theorem, 

lim E^iXgiS!")] = E^iXgiS-)]. 

m— f oo 

The rest of the proof is to verify E^ [Xg{S-)] — i which follows by a reverse 
use of Girsanov's transformation in Theorem 12^ and we leave it to the reader. 
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Figure 1: Approximation of the option price for various time steps N. The left 
panel takes T g [0, 1]; the right panel takes T G [1, 2]. 



— splitting ap proximal i or 



Figure 2: Comparison of the explicit price and the approximated price via 
splitting when 0, 5*2 = 1400, and N = 11. 
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Figure 4: Vulnerable option price against the writer's asset value . 




Figure 5: Impact of risk aversion and correlation. The left panel gives the option 
price against maturity for various risk aversion parameters 7. The right panel 
gives the price against various correlation parameters. We set /xi =0.1 and 
= 0.06 to satisfy the parameter restriction of Proposition 13.51 
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